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ABSTRACT 

Three years of optical monitoring of the low-mass X-ray binary (LMXB) 4U 1957-1-11 is 
presented. The source was observed in V, R and i-bands using the Faulkes Telescopes 
North and South. The light curve is dominated by long-term variations which are 
correlated (at the > Scr level) with the soft X-ray flux from the All-Sky Monitor on 
board the Rossi X-ray Timing Explorer. The variations span one magnitude in all 
three filters. We find no evidence for periodicities in our light curves, contrary to a 
previous short-timescale optical study in which the flux varied on a 9.3-hour sinusoidal 
period by a smaller amplitude. The optical spectral energy distribution is blue and 
typical of LMXBs in outburst, as is the power law index of the correlation /3 = 0.5, 
where Fi,_opt oc F^. The discrete cross-correlation function reveals a peak at an X-ray 
lag of 2-14 days, which could be the viscous timescale. However, adopting the least 
squares method we find the strongest correlation at a lag of ± 4 days, consistent 
with X-ray reprocessing on the surface of the disc. We therefore constrain the optical 
lag behind X-ray to be between -14 and -1-4 days. In addition, we use the optical- 
X-ray luminosity diagram for LMXBs as a diagnostic tool to constrain the nature of 
the compact object in 4U 1957+11, since black hole and neutron star sources reside 
in different regions of this diagram. It is found that if the system contains a black 
hole (as is the currently favoured hypothesis), its distance must exceed ~ 20 kpc for 
the optical and X-ray luminosities to be consistent with other soft state black hole 
systems. For distances < 20 kpc, the data lie in a region of the diagram populated 
only by neutron star sources (black hole systems are ten times optically brighter for 
this X-ray luminosity). 4U 1957+11 is unique: it is either the only black hole LMXB to 
exist in an apparent persistent soft state, or it is a neutron star LMXB which behaves 
like a black hole. 

Key words: accretion, accretion discs, black hole physics. X-rays: binaries, stars: 
individual: 4U 1957+11 



1 INTRODUCTION 

4U 1957+11 (V1408 Aql) is a persistently active (~ 
1-10x10"^" erg cm-2 s'^ 1.5-12 keV) low-mass X- 
ray binary (LMXB). Its X-ray s pectrum has always 
appeare d soft since its di scovery jGiacconi et al.l 1 19741 : 
[Wiinan ds. Miller fc van der Klis 200^. These two proper- 
ties make this source unique; most LMXBs are transient 
or regularly perform transitions to harder X-ray states 
(for more on X-ray state s see e .g. Hasinger fc van der Klis 
1 19891 : iHoman fc Bellonil l2005l : iMcClintock' fc RemillardI 



l2006l : iFender. Homan fc Bellonil |2009'). The X-ray spec- 
trum of 4U 1957+11 can be described by a blackbody 
plus a power law, wit h fluxes which vary on tim e scales 
of hours to months ([Y aqoob, Ebisawa fc Mitsuda' Jl993l : 
ISingh. Apparao fc KrafS , 1994 ; Nowak fc Wilm s 199^). An 
apparent 117 day X-ray period ( Nowak fc Wilms|[l999l ) pos- 



sibly due to a precessing accretion disc viewed almost edge- 
on was later shown to be much more complex; the long-term 
X-ray flux variations ar e likely to be driven by changes in 
the m ass accretion rate l| Yaqoob et al.|[l993l : I Wiinands et al.l 
l2002h . 
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No estimates of the mass function (e.g. ICasaresI |2007| ) 
have been established because the system has never faded 
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Figure 1. Finding chart for 4U 1957+11, indicating the two stars 
used for flux calibration. The image is a 100 sec exposure in i- 
band taken with the EMOl camera on FTN under good seeing 
conditions. The field-of-view is 60 X 60 arcsec. The optical coun- 
terpart of 4U 1957+11 is marked in the centre. The magnitudes 
of stars 1 and 2 are given in Table 2. North is to t he top and east 
is to the left. For a larger field, see iMareon et al] (1978). 



to quiescence since its discovery, and type I X-ray bursts 
typical of neutron stars have not been detected - as such it 
is uncertain whether the compact object in the 4U 1957+11 
system is a black h ole or a neutron star. It was argued 
jYaaoob et al.lll993l ) that the X-ray luminosity is too low, 
and the inner disc radius too large for a candidate black 
hole X-ray binary (BHXB), but this may not be the case, 
and several X-ray spec tral and timing properties are sim- 
ilar t o other BHXBs (jWiinands et al.l l2002l : iNowak et aD 
I2OO8I ). The X-ray spectrum is typical of soft state BHXBs, 
and its behaviour sha res many similarities with LMC X-3 
(jWiinands et al.ll2002l . and references therein), which con- 
tains a > 5.8 Mq black hole and is also in a persistent 
soft state (but unlike 4U 1957+11, L MC X-3 occasionall y 
performs transitions to the hard state; IWilms et al1l200lf ). 
Moreover, the luminosity is higher and consistent with soft 
state BHXBs if the distanc e to the source is further than the 
often assumed 7 kpc (e.g. iMarscon. Thorstensen fc Bowveil 
1 19781 ). The X-ray spectrum also appears slightly harder at 
higher fluxes, which is not usually the case for neutron star 
X-ray binarie s (NSXBs; Wiinands ct al. 20o2). 

Recently, iNowak et all (|2008l ) fitted high-resolution X- 
ray spectra of 4U 1957+11 with an inclined disc blackbody 
with an interstellar absorption of uh = (1.5 ± 0.5) x 10^^ 
cm~^. The total absorption through the Galaxy in the direc- 
tion of the source is 1.2 x 10^^ cm~^ (Kalberla ot al. 2005), 
so the measured absorption implies a large distance - on 
the far side of t he Galax y or possibly i n the halo (see also 
lYao et al.ll200i ). Indeed. [Nowak et all (|2008l ) use Chandra, 
XMM-Newton and Rossi X-ray Timing Explorer (RXTE) 
data to infer a rapidly spinning black hole and a distance 
of ~ 10-22 kpc. However, the authors do not explore the 
possibility of a neutron star using their data; instead they 
use the spectral and timing arguements above to favour a 
black hole. 



An optical counterp art to 4U 1957+11 was discovered 
by iMargon et al.l (|l978l ). which was found t o vary sinu- 
soidally on an orbital period of 9.33 hours l|Thorstensenl 
[l983). However, a later photometric study revealed a 
more complex light curve of larger amplitude variability 
iHakala. MuhU fc Dub"u3 Il999l ). In this latter study the 
authors establish a colour-dependency in the amplitude; 
the optical (UBVRI) spectra l energy distr ibution (SED) is 
slightly bluer at lower fluxes. iHakala et al.l (|1999. ) modelled 
the light curve assuming it modulates on the orbital pe- 
riod, however the data cover just one 9-hour orbit. Since 
the optical behaviour is in contrast with th e original sinu- 
soidal modulation found (|Thorstenseij 1987 ), the var i abilit; 
they measure may not be periodic. IShahbaz et al.l (|l99' 
took a broad-band (4000-10, OOOA) optical spectrum of 4U 
1957+11, which displayed a mostly featureless continuum, 
with weak Hq, H/3 and He II emission lines. No absorption 
lines from the companion star were detected. 

Here, we present the results of an optical monitor- 
ing campaign of 4U 1957+11 spanning 37 months during 
2006-2009. Concurrent X-ray monitoring reveals a positive 
optical-X-ray correlation, and helps to constrain the origin 
of the optical variability and the nature of the compact ob- 
ject. The observations, reduction and flux calibration are de- 
scribed in Section 2. We analyze the results and give a com- 
prehensive discussion, making comparisons to other sources, 
in Section 3. In Section 4 we summarize the results and con- 
clusions. 



2 OBSERVATIONS 
2.1 Optical photometry 

Observations of the field of 4U 1957+11 were taken with 
the Faulkes Telescope North situated at Haleakala on 
Maui, USA (FTN) and Faulkes Telescope South situated 
at Siding Spring in Australia (FTS). Both are 2- metre 
robotic telescopes optimized for research and ed ucation 
(|Roberts. Roche fc Rossll2008l : iLewis fc Rochell2009D . Imag- 
ing was obtained in Bessell V, Bessell R and Sloan Digital 
Sky Survey (SDSS) i filters between April 2006 and May 
200 9, as part of a mo nitoring campaign of ~ 30 LMXBs 
(see lLewis et al.ll2008l ). The log of observations is given in 
Table 1. Various cameras were used (column 3 of Table 1); 
all cameras have 2048 x 2048 pixels binned 2x2 prior to 
readout into effectively 1024 x 1024 pixels. For all cameras, 
the field of view is 4.7 x 4.7 arcmin and the pixel scale is 
0.278 arcsec pixel"^. Automatic pipelines de-bias and flat- 
field the Faulkes science images using calibration files from 
the beginning and the end of each night. 

The seeing ranged from 0.7 to > 3 arcsec. We discarded 
images that were taken when the seeing was greater than 
six pixels (1.67 arcsec), because the measured magnitudes of 
field stars of similar magnitude to 4U 1957+11 became inac- 
curate under those poor seeing conditions. We also discarded 
images with poor signal-to-noise ratio (S/N; mostly due to 
thin cloud) and images with tracking problems (where the 
stars are non-spherical and stretch more than a few pixels 
during the exposure), an excess of hot pixels, or a back- 
ground gradient across the image. 144 images were kept, 
which were taken between April 2006 and May 2009 (1138 
days). The discarded images are not listed in Table 1. 
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Table 1. Faulkcs Telescopes observations of 4U 1957+11 used in this work. 



Date MJD Camera Airmass Exposure times (sec) 

V R i 



2006-04-08 


53833.6 


DillCam 


1.26 




120 




120 


120 




2006-05-10 


53865.6 


DillCam 


1.08 




200 




200 


200 




2006-05-21 


53876.6 


DillCam 


1.04 




200 




200 


200 




2006-05-27 


53882.6 


DillCam 


1.02 




200 




200 


200 




2006-06-07 


53893.6 


DillCam 


1.01 




200 




- 


- 




2006-11-14 


54053.2 


HawkCam 


1.24 




200 




_ 


200 




2007-04-22 


54212.5 


HawkCaml 


1.56 




100 




- 






2007-06-16 


54267.5 


HawkCaml 


1.09 




100 




100 


100 




2007-07-10 


54291.4 


HawkCaml 


1.20 




100 




100 


100 




2007-07-12 


54293.5 


HawkCaml 


1.01 




- 




- 


100 




2007-08-06 


54318.3 


HawkCaml 


1.12 




100 




100 


100 




2007-08-10 


54322.3 


HawkCaml 


1.09 




100 




100 


100 




2007-08-19 


54331.4 


HawkCaml 


1.03 




100 




100 


100 




2007-09-02 


54345.3 


HawkCaml 


1.02 




100 




100 


100 




2007-10-05 


54378.3 


HawkCaml 


1.08 




100 




100 


100 




2007-11-13 


54417.3 


HawkCaml 


2.42 




100 




100 


100 




2008-06-13 


54630.4 


EMOl 


1.70 




100 




100 


100 




2008-06-16 


54633.5 


EMOl 


1.02 




80 




110, 100 


2 X 


100 sec 


2008-06-20 


54637.6 


EAOl 


1.12 




100 




100 


100 




2008-06-24 


54641.5 


EMOl 


1.05 




- 




2 X 100 sec 


2 X 


100 sec 


2008-06-25 


54642.5 


EMOl 


1.01 




- 




100 


2 X 


100 sec 


2008-06-30 


54647.6 


EA02 (FTb) 


1.20- 


-1.42 








18 X 100 sec 


2008-07-11 


54658.5 


EMOl 


1.01- 


-1.04 


2 X 


100 sec 


2 X 100 sec 


16 X 100 sec 


2008-07-14 


54661.5 


EMOl 


1.04 




100 




100 


100 




2008-07-14 


54661.6 


EA02 (FTS) 


1.37 




2 X 


100 sec 


2 X 100 sec 


4 X 


100 sec 


2008-07-16 


54663.5 


EMOl 


1.19 




100 




100 






2008-07-16 


54663.6 


EA02 (FTS) 


1.37- 


-1.42 


2 X 


100 sec 


100 


3 X 


100 sec 


2008-07-30 


54677.4 


EA02 (FTS) 


2.07 




200 




200 


200 




2008-08-09 


54687.3 


EMOl 


1.05 




100 




100 


100 




2008-08-31 


54709.5 


EMOl 


2.06 




100 




100 


100 




2008-09-10 


54719.5 


EMOl 


2.16- 


-2.41 


2 X 


100 sec 


100 


2 X 


100 sec 


2008-09-15 


54724.5 


EMOl 


2.13 








100 


100 




2008-09-25 


54734.4 


EMOl 


1.52 




100 




100 


100 




2008-11-04 


54774.2 


EMOl 


1.18 




100 




100 


100 




2009-02-22 


54884.6 


EMOl 


2.50 




100 




100 


100 




2009-05-17 


54971.5 


EMOl 


1.12 




100 




100 


100 





MJD = Modified Julian Day. The EA02 camera is on Faulkes Telescope South and all the other cameras are on Faulkes Telescope 

North. Data from 34 dates were used, totalling 144 images. 

We performed aperture photometry of 4U 1957-1-11 and 
two field stars using the package APPHOT in IRAF . Point- 
spread-function fitting was not used because minor tracking 
errors on some dates resulted in the stars being non-circular. 
A fixed aperture radius of 6 pixels was adopted (accurate 
centres to the stars are fitted), and a background annulus of 
10-20 pixels was used. 

Flux calibration was achieved using the standard star 
PG 1525-071 {V = 15.053; R = 15.141; / = 15.218) from 
the list of Landolt photometric standards (lLandoltlll992l ). 
which is observed on FTN regularly in V , R and i filters. 
On 2006-04-08, 2007-03-22 and 2009-02-22 the conditions 
were photometric and both PG 1525-071 and 4U 1957-f 11 
were observed by FTN. The SDSS i magnitude of PG 1525- 
071 was calculated (i = 15.53) adopting th e prescription de- 
scribed in Table 3 of IJordi. Grebel fc Amm on (200^) from 
the known R and / magnitudes. Taking into account the 
airmass differences between standard star image and tar- 
get image on each of the three nights, we calculated the V , 
R and i magnitudes of two field stars in the target images. 



Table 2. Measured magnitudes of two nearby field stars (Fig. 1) 
used for flux calibration. 

Star V R i 

1 16.47 ± 0.08 16.01 ± 0.08 15.85 ± 0.04 

2 15.94 ± 0.08 15.35 ± 0.08 15.10 ± 0.03 



The stars are marked 1 and 2 in Fig. 1 and their magnitudes 
are listed in Table 2. Errors were estimated from the range 
in values from the three dates. Finally, to convert appar- 
ent magnitudes to intrinsic, de-reddened flux densities, we 
adopt Ay = 0.93 for the in terstellar extinction towards the 
source (|Margon et al.l 19781) and u se the extinction curve of 



ICardeUi. Clavton fc Mathig (|l989l ) 



2.2 X-ray monitoring 

The All-Sky Monitor (ASM) on board the RXTE satellite 
monitors the X-ray sky several times per day between 1.5 
and 12 keV. 4U 1957-1-11 has been detected by the ASM 
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Table 3. Long and short-term flux and variability properties. 



Waveband 


Time span 


Mean counts 


1 (T 


Mean relative 


Mean unabs- 


1 cr 


Fractional 




used 


or magnitude 


range 


error 


orbed flux 


range 


rms 








4U 1957+11 long-term variability 






1.5-12 


13.5 


2.33 ASM cps 


± 1.03 


± 0.42 


6.9x10-1° 


± 3.1x10-10 


44 ± 18 % 


keV 


years 




ASM cps 


ASM cps 


ergs cm"'^ s" ^ 


ergs cm" ^ s^^ 




V 


3.1 years 


18.93 


± 0.23 


± 0.04 


0.241 mjy 


± 0.050 mJy 


21 ± 4 % 


R 


3.1 years 


18.77 


± 0.22 


± 0.04 


0.170 mJy 


± 0.034 mJy 


20 ± 3 % 


i 


3.1 years 


18.88 


± 0.21 


± 0.04 


0.176 mJy 


± 0.033 mJy 


19 ± 4 % 








Field star 3 long-term variability 






V 


3.1 years 


20.13 


± 0.11 


± 0.10 


unknown extinction 


< 19 % 


R 


3.1 years 


19.04 


± 0.04 


± 0.04 






< 7 % 


i 


3.1 years 


18.42 


± 0.04 


± 0.03 






< 7 % 



4U 1957-t-ll short-term variability 
2-60 keV 2.0 years (From Wiinands et al. 2002) 1-2 % 

i 0.9 hrs" 19.18 ± 0.08 ± 0.07 0.133 mJy ± 0.010 mJy 7.3 ± 6.5 % (1.1 a) 

i 1.4 hrs*" 18.92 ± 0.05 ± 0.02 0.169 mJy ± 0.008 mJy 4.9 ± 2.2 % (2.2 a) 



Field star 3 short-term variability 
i 0.9 hrs" 18.40 ± 0.05 ± 0.04 unknown extinction 4.6 ± 3.4 % (1.4 cr) 

i 1.4 hrs* 18.41 ± 0.02 ± 0.02 " < 3.0 % 

The 2-60 keV data from lWiinands et al.l l l2002h are RXTE Proportional Counter Array data from 1997-11-26 to 1999-12-02 but the 
variability timescales tested are 0.1-1000 seconds. All optical data are from the Faulkes Telescopes observations, "data taken on 

2008-06-30; ''data taken on 2008-07-11. 



since the start of its operation in 1996. The X-ray source 
has remained persistently active over the full 13.5 years: the 
mean daily count rate is 2.3 counts sec^^ and the rms is 1.0 
counts sec"!. This range corresponds to an unabsorbed flux 
of (6.9 ± 3.1) xlO-i° ergs cm'^ s"^ (1.5-12 keV) assuming 
a blackbody at 1 keV (which is typical of an X-ray transient 
in the soft state) an d an absorption of uh = 1.5 x 10^^ cm~'^ 
l|Nowak et al.|[200^ . 86 per cent of ASM 1-day averages are 
3 a detections, so it may be that on some dates the fluxes 
are significantly lower than the above range. We wish to 
cross-correlate optical and X-ray fluxes, so in order not to 
favour the brighter epochs, we did not discard the non-3 a 
detections. Instead, we binned the ASM data into 3-day av- 
erages, 5-day averages and so on to achieve higher S/N. This 
has the effect of smoothing the fight curve; daily variability 
information is lost (for the analysis in Section 3 we study 
the strength of optical-X-ray correlations as a function of 
the number of ASM days used). 



relative errors are smaller. 4U 1957-1-11 is variable, spanning 
a range of almost one magnitude in all three filters. 

From these light curves we derive the long-term X-ray 
and optical variability properties (Table 3) . Over these long 
timescales (years) the fractional root mean square (rms) of 
the X-ray variability is ^ 44 per cent and in optical it is 
~ 20 per cent (in all three filters) . We may expect a smaller 
optical variability amplitude if the optical light arises from 
reprocessed X-ray photons on the surface of the accretion 
disc, as is often seen in LMXBs in outburst. The optical 
variability is real; the relative errors in the optical photom- 
etry (~ 0.04 mag) are much smaller than the measured rms 
of the variability (~ 0.2 mag). The long-term optical frac- 
tional rms of the field star is < 7 per cent; the errors in the 
magnitudes due to the photometry are similar to the ap- 
parent variability, so the field star is consistent with being 
non-variable. From Fig. 2 it appears the X-ray and opti- 
cal fiuxes could be correlated; we explore this possibility in 
Section 3.5. 



3 ANALYSIS AND DISCUSSION 
3.1 Long-term variability 

The optical and X-ray light curves are presented in Fig. 2. It 
is clear from the ASM light curves (upper two panels) that 
the X-ray fiux is dominated by long-term variations on tens 
to hundreds of days, but short-term variability on timescales 
of days or less is also present. Binning the ASM data into 7- 
day averages has the effect of smoothing the low-S/N 1-day 
averages light curve and improves the clarity of structure on 
tens-to-hundreds of days. In the lower three panels of Fig. 
2, the i, R and V apparent (absorbed) magnitudes of 4U 
1957-f 11 and a field star of similar magnitude are shown. 
The magnitude errors in the figure are the total errors; the 



3.2 Short-term variability 

IWiinands eraD l|2002l ) showed that the X-ray variability on 
0.1-1000 second timescales is low, with amplitudes between 
one and two per cent. We tested for optical variability on 
timescales of minutes on two dates in 2008 (Fig. 3) . We chose 
the i-band filter because the S/N is usually higher in i than 
in R and V for a given exposure time. On 2008-06-30, 18 
observations were made with FTS, with a time resolution of 
~ 170 seconds. The source was faint (i ~ 19.2) and the S/N 
relatively low (likely due to a high airmass of up to 1.42). 
As a result, the relative magnitude errors were similar to or 
greater than the variability amplitude on this timescale, for 
both 4U 1957+11 and field star 3. On 2008-07-11 the source 
was brighter (i ~ 18.9), the conditions were favourable and 
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53800 53900 54000 54100 54200 54300 54400 54500 54600 54700 54800 54900 



MJD 

Figure 2. X-ray (1.5-12 kcV) and optical long-term light curves of 4U 1957+11. In the lower three panels, the magnitudes of 4U 1957-1-11 
and a local field star are shown. The errors represent the total systematic and photometric magnitude errors. 



6 D. M. Russell et al. 



18 p 
18.2 - 
18.4 - 
18.6 - 
18.8 - 

19 - 
19.2 - 
19.4 - 

14 

18 p 
18.2 - 
18.4 - 
18.6 - 
18.8 - 

19 - 
19.2 - 
19.4 - 



4U 1957+11 
Field star 3 



^ >^ ^ t >^ 



+ + 



14.2 14.4 14.6 14.8 
Hours since 2008-06-30 00:00 UTC 

" ' ' ' 4U 1957-Hll ^-^ 
Field star 3 — ^ 



15 



10.4 10.6 10.8 11 11.2 11.4 11.6 11.8 12 
Hours since 2008-07-1 1 00:00 UTC 

Figure 3. j-band light curves with time resolutions of hundreds 
of seconds, on two dates in 2008. 



the airmass very low. 15 observations were made with FTN 
within 36 minutes (the time resolution was ~ 140 seconds), 
and one 49 minutes before. It appe ars the sourc e brigh tened 
during that time by ~ 0.1 mag. iMotch et aLl l| 19851 ') also 
found 0.1 mag variations of the optical counterpart in one 
hour. 

The properties of the optical and X-ray (from 
IWiinands et al.ll2002l . ten years previously) short-term vari- 
ability are also given in Table 3. The i-band variability on 
2008-07-11 has a fractional rms of ~ 3-7 per cent. This 
is greater than the X-ray variability measured on these 
timescales, 1-2 per cent, however the i-band variability is 
significant at only the 2.2 a level. We do not detect variabil- 
ity from the field star, which is brighter than 4U 1957+11 
in i-band. 



3.3 Periodic variability? 

We carried out a more intensive monitoring campaign of 4U 
1957+11 in June - September 2008, in order to search for 
periodic variability on the claimed orbital peri od of 9.329 
± 0.011 hours established in 1985 (ThorstcnsenI 119871 ). We 
define zero phase at MJD 54646.0, at the centre of our cam- 
paign. This does not correspond to a known physical phase 
orientation of the binary. The above error of the period esti- 
mate propagates to an uncertainty of ~ 100 seconds per day. 
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Faulkes Telescope (2008): i-band 
Faulkes Telescope (2008): V-band 
Thorstensen (1987): V-band fit 
Thorstensen (1987): all V-band data 




0.5 1 1.5 2 

Orbital phase (P = 9.329 +- 0.01 1 hours) 

Figure 4. Optical magnitude versus orbital phase (plotted twice 
for clarity), assuming the ephemeris of lThorstensenI (1987). Phase 
zero is arbitrary and does not correspond to a linown orientation 
of the binary. The V and i magnitudes from June 2008 to Septem- 
ber 2008 were used (outside this range the error of the orbital 
phase is too large). 



which results in an uncertainty in the relative phase which 
increases with distance in time from MJD 54646.0. 

The folded light curve is presented in Fig. 4. We find 
no evidence for a modulation on this period, and the range 
of magnitude s in 20 08 greatly exceeds that found 23 years 
previously bv lThorste nsen ( 1 98?! ). Instead, the range of mag- 
nitudes we m e asure is consi stent with those seen in 1996 by 
iHakala et all l|l999D . Since iThorstenseiil (|l987l ) discovered 
the periodicity from observations spanning 28 periods (11 
nights), it is likely that the period was real at that time. 
The behaviour of the optical v ariability has c l early changed 
since 198 5, as pointed out bv iHakala et al.l (|l999l ). In the 
paper of IHakala et akl their observations covered just one 
9.3-hour period, and it was assumed that the variations they 
witnessed were associated with this period. We suggest that 
the variability may not have been periodic at that time (see 
also Sections 3.5). 

To test for the existence of other periodic signals, we 
performed timing analysis on our optical light curves us- 
ing the STARLINK package PERIOD. No significant intrin- 
sic signals were found in the i, R ot V light c urves using 
either the Lomb-Scargle (|Press fc Rvbickil 19891. and refer- 
ences therein) or epoch-folding (e.g. Horne. Wade &: Szkodvl 
1980) techniques, which are both appropriate for unevenly 
sampled data. For the i-band, in which we have short-term 
variability observations, we are able to search a frequency 
range of ~ 10~*-10~^ Hz (trial periods between ~ 100 sec- 
onds and 1 year). The strongest (but not highly significant) 
peaks are at periods of 24.4 hours, 49.8 days and 427.4 days. 
The Faulkes telescopes often observed the target around the 
same time each day during the monitoring, which probably 
causes the apparent 24.4-hour period. By folding the light 
curve on 49.8 days and 427.4 days, we see that these periods 
are dominated by the sampling of the two short-term vari- 
ability studies, so they are also sampling aliasing artefacts. 
We find no evidence for optical periods of hundreds of days 
like t hose uncovered at X -ray energies (jNowak fc Wilms! 
I1999I : IWiinands et al.ll20()3 ). Follow-up dedicated monitor- 
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Figure 5. Optical (de-reddened) versus X-ray fluxes centred on the same day. The solid lines are the best power law fits; the power law 
index /3 is shown, where -F,y,oPT oc F^. On the right panel we use RXTE ASM 7-day means, which result in the strongest correlations 
(see text and Fig. 6). 



ing campaigns covering various timescales would be benefi- 
cial for more stringent period searcfies. 

3.4 Optical spectral energy distribution 

The average optical {V to i) spectral index over tfie three 
years of monitoring was a = -1-1.0 (calculated from the mean 
de-reddened flux densities in Table 3), where (x v"' . From 
all dates in which both V and i observations were made, the 
1 a range of spectral indices observed was a — -1-0.92 ± 0.25. 
Since the extinction towards the source is known to be low 
(we adopt Ay ~ 0.93; see Section 2.1), our estimates of a 
are not likely to suffer large inaccuracies from an uncertain 
extinction. Systematic errors in the flux calibration could 
however produce errors in a up to 0.33. 

This optical spectral energy distri bution ( SEP) 

is typical of a LMX B in outburst l|Hvnesl l2005l : 

iRussell. Fender fc Jonkeil |2007| ). It is not consistent with 
the canonical a = i multi-temperature thermal disc spec- 
trum, nor the Rayleigh-Jea, ns tail of a blackbody a = 2.0 
l|Frank. King fc Rainel l20o3 ) . A redder SED would be ex- 
pected if the emission were synchrotron in nature, as is seen 
in the optical/infr ared (IR) SEDs of some LMXB jets (e.g. 
IChatv et al.ll2003h . With a ^ we are likely sampling a 
region of the disc blackbody between the peak, which resides 
in the ultraviolet (U V), and the Rayleigh- Jeans tail in the 
IR (e.g. lHvnesllioOsI ). However, it is uncertain whether this 



is the Rayleigh- Jeans tail of the multi-temperature viscous 
disc (in which case a — ^ in the UV) or the irradiated disc 
(where the disc surface reprocesses th e X-ray photons to 
lower energies; here a = — 1 in the UV: lFrank et al.ll2002l ). 

We would expect the optical SED to steepen (values 
of a to increase) with disc temperature, for a disc of con- 
stant area. This was recently shown to be the case for 
th e optical-IR colours of a number of outbursts of Aql X- 
1 (|Maitra fc Bailvnl 1200^ ) and successfully modelled as a 
single-temperature heated blackbody. The V-i colour of 4U 
1957+11 is indeed bluer when it is brighter; in a forthcom- 
ing work we will apply the same technique to our data of 
4U 1957-1-11 and other LMXBs (Russ ell, Maitra et al. i n 
preparation; for preliminary results see lRussell et al.]|2008l ). 



3.5 Optical— X-ray correlations 

If the accretion disc (viscously heated or irradiated) is 
producing the optical emission, as implied by the SEDs, 
we expect the optical and X-ray hght curves to be cor- 
related. Changes in the mass accretion rate (and possi- 
bly a precessing disc) produce the long-term X-ra y varia- 



bly a precessing disc ) produce tne long-term A-ra y varia- 
tions ((Nowak &: Wilmj|l999l: IWiinands et alll2002l). It has 



been shown ( van Paradiis fc McClintocklll994l :l iRussell et al.l 
I2OO6I ) that optical-X-ray correlations exist for LMXBs over 
orders of magnitude in luminosity. The power law index of 
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the correlations depend on the process that produces the 
optical emission. 

For 4U 1957+11, we take the mean magnitude from 
each day (0-24 UT) and plot the optical flux density against 
the X-ray 1.5-12 keV count rate centred on the same day, 
in Fig. 5. A positive optical-X-ray correlation exists which 
approximates a power law. The scatter in the correlation is 
smaller if we use 5 or 7-day ASM average count rates in- 
stead of 1-day averages (shown for i-band in the left panels 
of the figure). This is probably because the S/N of the 1-day 
averages is low, and the X-ray light curve is dominated by 
variations on timescales longer than one day. We do not in- 
clude the errors associated with the optical and X-ray fluxes 
in fitting the correlations. 

We investigate the strength of the correlation as a func- 
tion of the number of days of ASM data averaged, in Fig. 6. 
The residuals to the fit (measured in log F;^,opt) are squared 
and summed, so the minimum in this plot corresponds to the 
least squares, and the strongest Fqpt-Fx correlation. The 
best fit is obtained using 7-day ASM averages, for V and 
i-bands, and for the mean of the V, R and i-bands. The 
residuals are approximately twice as large if 1-day averages 
are used, confirming the low S/N of these data. If > 9-day 
averages are used, the residuals increase again. This is likely 
to be an indication of the timescale of the X-ray variabil- 
ity - structure in the intrinsic X-ray light curve is diluted, 
or smoothed by averaging data over timescales longer than 
the variability. We therefore use 7-day ASM averages for the 
correlations presented in the right panel of Fig. 5. 

To test the significance of the correlations we make use 
of the Spearman's Rank correlation. A correlation exists be- 
tween FopT and Fx at the 4.4 a, 4.3 a, 4.5 a and 4.0 a 
levels for V, R, i and the mean of V ,R,i-ha.nds, respectively. 
The significance of the correlation is highest when i-band is 
used, for which we have the most data points, and least when 
the mean of the three filters is used, for which we have the 
least number of data points. Since the variability timescale 
is on the order of days, these significances may be overes- 
timated by many data points separated by less than this 
timescale. We therefore calculate the Spearman's Rank cor- 
relation again for i-band, removing all data within 20 days 
of an other data point. Data on 18 dates were left, for which 
a Spearman's Rank correlation exists with a significance of 
3.1 a. 

The fits to the correlations have power law indices 
of l3 — 0.5-0.6 (where Fqpt oc F^), which is empiri- 
cally typical for both BH XBs and NSXBs ov er several or- 
ders of magnitude (Russe ll et all |2006| . l2007h . Within er- 
rors, /3 is no different if V, R, i-band data or the mean 
of the three are used. P « 0.5 is expected if X-ray re- 
processing on the surface of the disc produce s the op- 
tic al emi ssion (van Paradiis fc McClintockl Il994l : although 
see iDubus . Hamcurv & Lasota 2001 for a relation with an 
evolving /3 based on modelling light curve morphologies). 
The same relation is also expected if the viscously heated 
disc dominates, but only in the c ase of a NSXB; /3 is ex- 
pected to be shallower for a BHXB fauss ell et aLll2006 l). Op- 
tical synchrotron emission from a BHX B jet would produc e 
a slightly steeper power law, /3 « 0.7 (|Russell et al.l [20061 ): 
although the measured power law index for 4U 1957-1-11 is 
close to this value, the spectral index is inconsistent with a 
synchrotron origin (Section 3.4). 
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Figure 6. The strength (least squares) of the -Fqpt^^x corre- 
lation versus the number of days of averaged RXTE ASM data 
used (centred on the date of each optical observation). The ver- 
tical line indicates seven days of averaged X-ray data, for which 
the strongest -Fqpt^^x correlations exist. 

So far we have investigated optical-X-ray correlations in 
which the X-ray data are centred on the same day as the op- 
tical. It is also possible to test whether one of the light curves 
lags the other. The disc rete cross-correlation function (CCF; 
lEdelson fc Kroliklll988l ) is the standard, effective method for 
two unevenly sampled light curves. In Fig. 7 we plot the dis- 
crete CCF for optical and X-ray fluxes, using i-band (upper 
panel) and the mean of V, R and z-bands (centre panel). We 
take the mean optical magnitude per filter on each day in 
which two or more observations were made. 

A positive correlation is indeed seen, peaking at an X- 
ray lag of ~ 7 days. We estimated the error associated with 
this lag thus. We removed one third of the i-band data (10 
random data points from the 32) and performed the discrete 
CCF again, measuring the peak lag. 25 iterations of this 
were carried out using i-band and ASM 7-day averages (first 
trial) and 25 further iterations using i-band and ASM 3- 
day averages (second trial). The mean lag and 1 a error 
were 7 ± 1 days for the first trial and 8 ± 4 days for the 
second trial. In none of the 25 iterations was an optical lag 
measured in either trial (the X-ray lag was always greater 
than zero). Since seven days of averaged ASM fiuxes were 
used for each optical data point in the first trial, the formal 
error on the measured lag in that trial cannot be less than 
±4 days. Likewise for the second trial the error on the lag 
must be no less than ±2 days. We therefore measure X-ray 
lags of 7 ± 5 days from the first trial and 8 ± 6 from the 
second trial. Using the most conservative of these errors, we 
derive a measured X-ray lag of 2-14 days from the discrete 
CCFs. 

For X-ray reprocessing, we would expect an optical lag 
on li ght-travel timescales (seconds; see e.g. lO'Brien et al.l 
120021 ) ■ i.e. consistent with zero days lag. The optical pre- 
cedes the X-ray by several days, which is indicative of op- 
tical variations from the intrinsic viscous disc. X-ray delays 
of several days are expected if matter or waves of turbu- 
lence propagate inwards from the optically-emitting outer 
disc to the X-ray-emitting inner disc. This implies that for 
4U 1957-1-11, accretion rate perturbations in the disc take 2- 
14 days to propagate from the optically-emitting outer disc 
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Figure 7. The discrete CCF using j-band/1.5-12 keV (upper 
panel) and mean(y,i?,i)/1.5— 12 keV (centre panel). Lower panel: 
the strength of the fit to the optical-X-ray correlation, as mea- 
sured from the sum of the squares of the residuals of the fit, versus 
the optical delay (see text). The error bar at the bottom of the 
plot represents the 7-day range of X-ray data used. 



to the X-ray-emitting inner regions close to the compact 

object. 

iBrocksopp. Groot fc Wiin3 l|200ll ) find a similar dis- 
crete CCF for the BHXB LMC X-3, with an X-ray lag of 
5-10 days from long-term optical and X-ray light curves. 4U 
1957-1-11 and LMC X-3 share many properties - both are 
persistent with a predominantly soft X-ray spectrum. It was 



found that both the optical and X-ray long-term variations 
of L MC X-3 are driven by changes in the ma ss accretion 
rate (|Wilms et al.ll200ll : IBrocksopp et al]|200ll ). We arrive 
at the same conclusion for 4U 1957-1-11. It is also worth not- 
ing that the X-ray light curve of the NSXB 4U 0614-1-09 
also varies with the mass accretion rate (|van Straaten et al.l 
[2000) and the optical variability may be correl ated with the 
X-rav (iMa chin ct al .ll990l : IShahbaz et al.ll2008h . In addition, 
iHoman et al.i (|2005i) found an X-ray lag of 15-20 days com- 
pared to near-IR variations in the BHXB GX 339-4 during 
a soft state. These measured lags probe the disc viscosity; 
they may be a function of both the size and viscosity of the 
disc. 

A further test of optical or X-ray lags is to directly 
calculate the strength of the optical-X-ray correlation as a 
function of the lag. We fit the Fqpt-Fx correlation tak- 
ing X-ray data (7-day averages) centred at a series of days 
before and after each optical observation. We found that 
the correlation was strongest (the sum of the squares of the 
residuals was least) when X-ray data centred on the same 
day as the optical were used (independently for all three op- 
tical filters). The error associated with this lag is ±4 days 
because ASM 7-day averages were used. The strength of the 
FoPT-Fx correlation is plotted against the optical delay in 
the lower panel of Fig. 7. The correlation is not strong for 
an X-ray lag of > 4 days. The discrete CCF peaked at an 
X-ray lag of 2-14 days; this lag is not reproduced here using 
the least squares method. 

The main difference between the two methods (dis- 
crete CCF and least squares to the Fqft-Fx correlation) 
is that the discrete CCF is normalised by the variance of 
the light curves (the product of the standard deviations). If 
the variability characteristics are very different on different 
timescales (as would be expected if there were two distinct 
processes contributing to the variability), the shape of the 
resulting correlation function will be more sensitive to the 
correlation that is strongest on that particular timescale. If 
this is the case then the discrete CCF is sensitive to the vari- 
ability on timescales of days (such as the viscous timescale) 
and the least squares method is most sensitive to timescales 
no greater than ~ one day (light travel timescales are <C 
one day). In addition, the least squares method tests for 
a power law relation between Fqft and Fx- This method 
would therefore be ineffective if the two fluxes are correlated 
but not by a single power law. 

Alternatively, the two methods may simply propagate 
errors systematically in different ways, resulting in different 
measurements for the same lag. The two methods result in 
lags that indeed overlap (both are consistent with an X-ray 
lag of 2-4 days) and since the systematic errors using each 
method are unknown, we must take into account all the 
errors (especially since the S/N of the X-ray data is low). 
We therefore conclude that the optical lag is constrained to 
be between —14 and +4 days. 

Signatures of both phenomena - variability from the 
viscous disc and the X-ray heated disc - are also found 
in the opti cal emission of active galactic nuclei (AGN; see 
iBreedt et~al..,2009 . and references therein). Here, the larger 
size scales of the accretion disc and supermassive black hole 
result in reprocessing timescales of days instead of seconds 
and viscous timescales of months-years instead of days. Like 
LMXBs, X-ray reprocessing tends to be the dominating 
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Figure 8. X-ray luminosity versus monochromatic optical luminosity for all BHXBs and NSXBs collected {^sseiTeTaTyioolJSoO^ 
The neutron star 'Z-sources' are plotted as separate symbols to neutron star 'atoll' sources (for definitions see Hasingcr fc van der Klij 
ll98Sl) and millisecond X-ray pulsars. The V- and i-band data for 4U 1957+11 are overplotted assuming three distances encompassing the 
minimum and maximum likely values. In the left panel all luminosities are plotted, from quiescence to peak outburst; the high luminosity 
region is expanded in the right panel. 



source of optical variability in AGN, but years of monitor- 
ing of some AGN have uncovered long-timescale variations 
from mass ac cretion rate changes or thermal fluctuations in 
the disc (e.g. IShemmer et al.ll20o"3 : iMarshall. Rvle k. Milieil 
I2OO8I : lArevalo et al.ll2009l '). One difference appears to be that 
the amplitude of the very long-term optical variations in 
AGN exceeds that of the X-ray, whereas in LMXBs (at least 
for 4U 1957-f 11 and LMC X-3) the long-term X-ray varia- 
tions are of higher amplitude than the optical. 



3.6 A black hole or a neutron star system? 

The nature of the compa ct object in th e 4U 19 57-1-11 system 
is disputed (Section 1). I Russell et al.l l|2006l ) showed from 
data collected from 24 LMXBs that BHXBs and NSXBs 
reside in different regions of the optical-X-ray luminosity 
diagram (Figures 1 and 2 in their paper). It would therefore 
be fruitful to utilize this as a diagnostic diagram to constrain 
the BHXB/NSXB nature of 4U 1957-Hl. In Fig. 8 we plot 
the quasi-simultaneous X-ray luminosity and optical (BVRI- 
bands; we exclude all near-IR data) monochromatic lumi- 
nosity (flux densit y scaled to distance) for a ll BHXBs and 
NSXBs collected in lRussell erahl l|2006l , l2007h . The distance 
to 4U 1957+11 is unknown, so we overplot the optical-X- 



ray data of 4U 1957+11 adopting three trial distances. The 
minimum distance is ~ 7 kpc and from high-resolution X-ray 
spectra the source is consistent with a B HXB at 10-22 kpc 
l|Margon et al.lll978l : pNowak et al.ll2008l ). If the s ource is in 
the G alactic halo its distance could be ~ 40 kpc (|Yao et al.l 
l2008h . 

At a distance of 7 kpc, the source lies amongst the 
NSXBs in the diagram and has an X-ray luminosity of 
^036-37 gj.g g-i (2-10 keV). AU BHXBs are ~ one order 
of magnitude brighter in optical at this X-ray luminosity 
(and are hard state BHXBs not soft state BHXBs), so this 
would strongly favour a NSXB. At 22 kpc, the data still lie 
in a region of the diagram populated by the NSXBs, but 
almost overlap with some of the soft state BHXB data (we 
should not compare with the hard state BHXB data because 
4U 1957+11 is a soft state source; see the right panel of Fig. 
8 for an expanded view of this region of the diagram). If 
4U 1957+11 is in the halo at 40 kpc, the data are typical 
of a (fairly optically faint) soft state BHXB, or a 'Z-source' 
NSXB. At 40 kpc the source may be super Eddington; for a 
10 Mq black hole the bolometric luminosity observed would 
be twice the Eddington luminosity. We note that LMC X-3 
is optically brighter than 4U 1957+11 (Fig. 8; right panel). 



An optical-X-ray correlation in 4U 1957+11 11 



We propose three scenarios for the nature of 4U 
1957+11: 

• It is a BHXB on the far side of the Galaxy at d~ 20- 
25 kpc. This would be consistent with both the X-ray 
timin g/ spectral evidence for a BHXB (e.g. IWiinands et alj 
I2OO2I ) and (marginally) the optical/X-ray luminosities of a 
soft state BHXB (Fig. 8); 

• It is a halo BHXB at d > 25 kpc. The optical/X- 
ray luminosities would be typical of other soft state 
BHXBs in this case but BH XBs in the halo are unusual 
(IWhite fc van Paradiis|[l996l ). This scenario would have im- 
portant implications for models of super nova 'kicks' and 
BHXB formation Ce.g. iMirabel et al.ll200ll '): 

• The source is in the Galaxy, at d;* 7-20 kpc. It would 
then be peculiar: it would either be an anomalous NSXB 
with X-ray properties like no other NSXB, or the faintest 
soft state BHXB, with the lowest normalization in the 
optical-X-ray luminosity diagram. 



4 CONCLUSIONS 

We have monitored the optical flux of the persistent LMXB 
4U 1957-1-11 in three wavebands over a period of three years 
with the two 2-m Faulkes Telescopes. Like the soft X-ray 
light curve, the optical continuum is dominated by long- 
term variations on timescales of days-months, likely origi- 
nating in changes of the mass accretion rate. The amplitude 
is a factor ~ two in flux over these timescales; around half 
the amplitude of the X-ray variations. Contrary to previous 
shorter-timescale studies, no clear periodicities exist in the 
optical light curves. However we do not dispute the orbital 
period of 9.3 hours previously identified over a time span 
of ~ 10 days bv [rhorst enscn (1987). We suggest the X-ray 
flux (and so X-ray heating of the disc and star) during their 
observations may have been fairly stable (as is sometimes 
the case on those timescales; see Fig. 2) so that the low- 
amplitude sinusoidal variations from the heated face of the 
companion star were revealed (to our knowledge no X-ray 
data were acquired during that epoch). Under the assump- 
tion that the vari ability they se e over one 9-hour period in 
1996 was periodic. iHakala et al.l |l999) model their complex 
light curve as an evolving accretion disc, and infer a mass 
ratio q ~ 0.4-0.5 and an inclination angle of 70-75°. We 
remark that these estimates of the system parameters are 
invalid as the flux was probably not periodic at that time. 

The optical SED is blue (a « -1-1.0) and slightly redder 
at lower fluxes. This is indicative of an accretion disc (heated 
viscously or by irradiation) and rules out synchrotron emis- 
sion (which has an even redder SED) as the dominating pro- 
cess. We establish a positive correlation between the optical 
and RXTE ASM 1.5-12 keV X-ray fluxes, of F^.opt oc 
significant at up to the 4.5 a level (conservatively, > Scr) 
and evident from the three optical wavebands independently. 
The power law index of the correlation shows that the em- 
pirical correla tion found over whole outbursts for LMXBs 
jRussell et al . 2006) can hold over much smaller ranges in 
luminosity. 

Cross-correlating the optical and X-ray light curves 
with the least squares method reveals a lag of ± 4 days, 
consistent with X-ray reprocessing in which X-ray variations 



should precede the optical on light-travel timescales (sec- 
onds). However, the peak in the discrete CCF is at an X-ray 
lag of 2-14 days. This could be the viscous timescale for 
matter or waves of turbulence to travel from the optically- 
emitting outer disc to the X-ray-emitting inner disc. Com- 
bining these two methods we constrain the optical lag to 
be between -14 and +4 days. Optical and X-ray monitoring 
with more data and with higher S/N in X-ray is required to 
further constrain these lags. 

BHXBs and NSXBs lie in different regions of the 
optical-X-ray luminosity diagram. We therefore overplot the 
data of 4U 1957+11, adopting a range of distances to the 
source, in order to constrain the currently disputed nature 
of the compact object. At distances < 20 kpc the data lie 
in a region of the diagram populated only by neutron star 
sources. If it were a BHXB at this distance, the optical emis- 
sion would be one order of magnitude fainter than the other 
BHXBs, and it would be the faintest soft state BHXB. We 
therefore postulate that if its distance is < 20 kpc it is most 
likely a peculiar neutron star system. 

However, taking into account the BHXB-like X-ray 
properties of 4U 1957+11, the most likely scenario is that 
the source is a distant BHXB. Its optical and X-ray lumi- 
nosities are consistent with other soft state black hole sys- 
tems if its distance exceeds ~ 20 kpc. 4U 1957+11 may 
well be a Galactic halo BHXB, in analogy with the proto- 
typical halo BHXB, XTE J1118+480. This would imply a 
high space velocity and a large supernova 'kick' at its birth. 
Optical phase-resolved radial velocity studies during a faint 
period (when the star/disc ratio is highest) may be required 
to identify the compact object unambiguously. 

If 4U 1957+11 indeed contains a black hole, it is a valu- 
able addition to the very few persistently accreting black 
hole X-ray binaries known (e.g. LMC X-1; LMC X-3) and 
possibly only the second in our Galaxy, after Cygnus X-1. 
In addition, 4U 1957+11 would be the perfect candidate to 
test for the presence of a weak soft state radio jet. A steady, 
core radio jet has never been dete cted from a BHXB in the 
soft state (e.g. iFender et al.|[2009l ) but it is speculated that 
faint jets (either intrinsically weak or with very low radia- 
tive efficiency) may be produced by soft state BHXBs. The 
only radio observations of 4U 1957+11 in the literature pro- 
vide an upper limit of 2.1 mjy at 5 GHz l|Nelson &: Spencen 
Il988l ). X-ray binaries can be detected down to, for example 
50 /i Jy or fainter with the Very Large Array (iGallo et al.l 
I2OO6I ) and much fainter with future planned Square Kilome- 
ter Array Pathfinders. 
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